The metallicity dependence of WR wind models 



G. Grafener & W.-R. Hamann 

Department of Physics, University of Potsdam, Germany 

Abstract. With the advance of stellar atmosphere modelling during the last 
few years, large progress in the understanding of Wolf-Rayet (WR) mass loss has 
been achieved. In the present paper we review the most recent developments, 
including our own results from hydrodynamic non-LTE model atmospheres. In 
particular, we address the important question of the Z-dependence of WR mass 
loss. We demonstrate that models for radiatively driven winds imply a rather 
strong dependence on Z. Moreover, we point out the key role of the L/M-ratio 
for WR-type mass loss. 



1. WR-type stellar winds 

WR stars show exceptionally strong winds with mass loss rates of the order 
of the single-scattering limit or above. The observed wind efficiency numbers 
77 = Mvoo/iL+fc), which denote the ratio between the wind momentum and the 
momentum of the radiation field, are typically in the range of 1-5. Therefore, 
if WR-type winds are driven by radiation, photons must be used more than 
once, i.e., after absorption photons must either be re-distributed into different 
wavelength regimes or they must be scattered more than once by overlapping 
lines in the extended wind. Because of their high wind densities, WR stars de- 
velop extended atmospheres where the hydrostatic layers are completely hidden 
to the observer. In combination with their strong radiation fields, they develop 
extreme non-LTE conditions that require sophisticated modeling techniques. 

The key question concerning WR mass loss is why a star becomes a Wolf- 
Rayet star. The observed parameters of early- type O stars, namely their lumi- 
nosities and effective temperatures, are in principle comparable with late- type 
WN stars. O supergiants even tend to show a similar He- and N-enriched sur- 
face composition. Because of their higher mass loss rates, however, WR stars 
show dramatically different spectra with strong and broad wind emissions where 
O stars show tiny photospheric absorption lines. 

If the WR winds are radiation-driven like OB star winds, one would also 
expect a similar Z-dependence of their mass loss. In this case, however, the 
question must be addressed why WNL stars are so different from O stars. 



2. Monte-Carlo models 



The first wind models which were able to reproduce the high efficiency factors 
of WR winds wer e obtained by Monte-Carlo techniques ( Lucv h Abbottlll993l : 
ISpringmannlll994l ) . In this approach the path of single photon energy packets, 
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released at the wind base, is followed until they escape from the wind. On their 
way through the wind, the photons are scattered by Doppler-shifted spectral 
lines which are treated as infinitely narrow pure scattering opacities. For each 
scattering process the energy is determined which is transferred from the radi- 
ation field to the wind material or vice versa. Finally a 'global' solution for the 
mass loss rate is obtained by comparing the total mechanical wind energy which 
is gained from the radiation field per time interva l with the overall wind lumi- 
nosity M«^. In the approach of iLucv h Abbottl ( 1993h the line opacities are 
taken from a line list and the atomic populations and ion ization are calculated 
in a modified nebular approximation. ISpringmannl l| 19941 ). on the other hand, 
used a purely statistical approach for the strength and the distribution of the 
line opacities. 

The pilot studies bv lLucv fc Abbott! (jl993) and ISpringmannl jl994) showed 
that high wind momenta with r\ up to 10 can in principle be achieved by multiple 
line scattering if enough line opacities are present. Nevertheless, these models 
do not explain why the WR mass loss rates are so high. They rely on an adopted 
/3-type velocity structure with a prescribed terminal wind velocity Voo (note that 
the wind luminosity goes with v^). In particular, the models fail to reproduce 
the wind acce leration in the deeper at mospheric layers below a radius of ~ 2R*. 

Recently, IVink de Koterl ((20051 . see also this volume) inv estigated the Z- 
depen dence of WR mass loss by means of similar models as ILucv fc Abbottl 
( 1993). In their approach, however, the excitation temperatures which enter the 
calculation of the ionization structure a nd the atomic pop ulations, are extracted 
from a grid of model atmospheres (see IVink et al.lll9 99). For the example of 
a late-type WN star IVink k, de Koterl find a similar Z-dependence as for O- 
stars with M oc Z°' m which flattens at metallicities below 10~ Zq because the 
line driving is taken over by N and He. Note that the mass loss rates at this 
metallicity are so small (~ lO~ 7 M0/yr) that the star would not be identified as 
a WR-type from its spectral appearance. For a WCL model they find a relation 
with Z om which flattens already at W^Zq. 



3. Optically thick wind models 

A approach which is complementar y to the previous one is the critical-point 
analy sis for optically thick winds (e.g. lPistinner Eichlerll995l : lNugis &: Lamerd 
2002). The underlying assumption of this method is that, because of their high 
mass loss rates, WR atmospheres are so extended that the sonic point of the 
wind flow is located at large flux-mean optical depth t s . If r s is large enough 
the radiation transport can be treated in the diffusion limit. The expression for 
the radiative acceleration then simplifies to 

Orad = - / XvF v dv = XRossT— y" i 1 ) 

c J 4:7rr z c 

where XRoss denotes the Rosseland mean opacity which ca n easily be obtained 
e.g. from the OPAL opacity tables ( Iglesias fc Bogerslll99fi[ ). This expression for 



a ra d is not only easy to evaluate, it also changes the dynamic wind properties 
in relation to OB stars. The wind acceleration in a spherically expanding radia- 
tively driven flow is given by the contributions of the gravitational attraction, 
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Figure 1. Solution of Eq. JSJ) in the p-T plane. The sonic-point conditions 
for an optically thick wind, i.e. xr oss = Xcrit with outward increasing XRoss, 
are fulfilled on the solid parts of the curve between 40 and 100 kK, and above 
160 kK. The Rosseland opacities in this plot are taken from the OPAL opacity 
tables. Typical WC star parameters are assumed. 



the pressure gradient, and the radiative acceleration 

dv MG I dp 

"J- = o l-+Orad- (2) 

dr r z p dr 

When the velocity gradient is extracted from the pressure gradient, this equation 
becomes 

/ a 2 \ dv MG a 2 da 2 

[v — = =- + 2 — + a rad 3 

i v I dr r z r dr 

where a denotes the sonic speed. In the diffusion limit a rad does not depend on 
4^. Eq. © then has a critical point at the sonic radius r s where v = a. A finite 

value of 4 s can only be obtained if the right hand side of Eq. © is zero at this 
point. 

MG a 2 da 2 

= + 2 — + a rad 4 

r| r s dr s 

In contrast, for the thin winds of OB stars a ra( j depends on 4^ via Doppler 
shifts. The critical point of Eq. © is then located at a significantly higher speed 
(the so-call ed Abbott sp eed) which corresponds to the fast radiative-acoustic 
wave mode (lAbbottJ ll980b For optically thick winds Eq. @ serves as a limiting 
condition for the mass loss rate. For reasonable wind parameters the second and 
third term on the right-hand side of Eq. (J3J) become negligible so that 

-jT ~ o rad (r s ) = Xcrit -^r c ■ (5) 
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The Eddington limit, referring to the Rosseland mean opacity, is therefore 
crossed at the sonic point. The critical value Xcrit for the Rosseland mean opacity 
can be calculated directly from the given L/M ratio for a specific object. 

In Fig. ^ we plot the the solution of Eq. (JJjJ), i.e. the relation between density 
and temperature where XRoss(p,T) = Xcrit> for the example of a typical WC 
star. In the hydrostatic layers below r s , the radiative force must be lower than 
the Eddington value. XRoss therefore has to increase outward, with decreasing 
density. As shown in Fig.^this condition is fulfilled at the hot edges of the two Fe 
opacity peaks, at temperatures in the range of 40-100 kK and above 160 kK. The 
resulting mass loss rates on these parts of the curve are given by M = AtrR+pa, 
where the sonic speed a depends on T and the chemical composition of the wind 
material. 

To estimate the actual density and temp erature at the sonic poin t, however, 
a further constraint is needed. At this point Nu eis Sz Lamersl (I2002T) utilize an 
approximate relation between temperature and optical depth bv iLucvl ( 197lf ) 



T\r) = (V(r) + |w(r)) (6) 

with the modified optical depth r' and the dilution factor W which is close to 
unity in our case, r' is obtained from the assumption that the outer wind is 
driven by radiation. For this purpose it is necessary to take the density and 
velocity structure above the sonic point into account. Again, the results rely 
on the adopted velocity distribution v(r), where a /3-type velocity law with a 
terminal wind veloci ty in the observed range of 1000-3000 km/s is used. 

The analysis by (Nueis &: Lamersl reveals important insights how the mass 
loss of WR stars might be adjusted. They find that the observed WR mass 
loss rates are in agreement with the optically thick wind assumption, with two 
distinct sonic-point temperature regimes. The cool regime with typical T s in 
the range of 40-100 kK and r s = 3 — 10 corresponds to late-type WN stars, and 
the hot regime with T s ~ 160 kK and r s = 1 — 30 to early-type WN and WC 
stars. For the mass loss of early-type WN stars, Nugis (these proceedings) finds 
a Z-dependence with very steep exponents in the range 1-1.5, dependent on the 
stellar mass. When the mass loss goes down at low metallicities, the validity of 
the optically thick wind assumption is however questionable. 



4. PoWR hydrodynamic model atmospheres 

From Monte-Carlo models we have learned that the acceleration of the outer 
part of WR winds, namely their high wind performance numbers, can in prin- 
ciple be ex plained by multip le line-scattering. In addition, the critical-point 
analysis by iNugis &; Lamersl revealed that the mass loss rates are presumably 
adjusted by the radiative driving of Fe-peak opacities at large optical depth. 
As a consequence we expect two regimes of WR mass loss corresponding to the 
two Fe-opacity peaks at cool and hot temperatures. Moreover, in addition to a 
Z-dependence, we expect a strong dependence of M on the L/M ratio. 

Regardless of these important findings, fundamental questions still remain 
open. All previous models rely on pre-defined velocity distributions v(r). It is 
therefore not clear if the high mass loss rates can be maintained throughout the 
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whole wind by radiation pressure alone. In reality, the terminal wind velocity 
adjusts by an interplay between the acceleration in the outer parts of the wind 
and the limiting condition at the critical point. Moreover, the optically thick 
wind assumption, i.e., the question whether the diffusion limit is valid at the 
critical point or if the radiation field is affected by Doppler shifts, remains to be 
verified. Finally, spectral analyses and time resolved spectroscopy show direct 
evidence for a clumped wind structure, but the effects of clumping are not 
included in any previous stellar wind models. 

The new Potsdam Wolf-Rayet (PoWR) hydrodynamic model atmospheres 
fseelGrafener h Hamannl2005l : lHamann & Grafener 2n0.i lKoeste'rke" et al . 2002; 
IGrafener et al.ll2002D address exactly the questions explained above. By a com- 
bination of lin e-blanketed non-LTE model atmospheres with the equations of hy- 
drodynamics, |&lfen£^^^Hamann| ( 2005h obtained the first fully self-consistent 



models for WR winds. In these models M, v(r), T(r), and the full set of non-LTE 
population numbers are computed consistently with the radiation field in the co- 
moving frame (CMF), i.e., the radiation transport, the equations of statistical 
equilibrium, the energy equation, and the equation of motion are simultaneously 
solved. Through the exact solution of the radiation transport no simplifying as- 
sumptions concerning the radiative acceleration are necessary. a ra d is obtained 
from integrating the product of opacity and flux over frequency (see Eq.^) , which 
are both calculated on a fine frequency grid in the co-moving frame of reference. 
Moreover, clumping is taken into account in the limit of optically thin clumps. 
The resultant models describe the conditions in WR atmospheres in a realistic 
way, and provide synthetic spectra, i.e., they allow a direct comparison with ob- 
servational material. Nevertheless, simplifying assumptions are still necessary. 
These are especially the assumption of a constant Doppler broadening velocity 
throughout the atmosphere, and the omission of opacities, partly due to the 
neglect of trace elements like Ne, Ar, S, or P and partly due to incompleteness 
of the available data. 

In the following we present PoWR models for both expected temperature 
regimes. 'Hot' models with core temperatures of T* = 140 kK for early-type 
WC stars on the He-main sequence, and 'cool' models for late-type WN stars 
in the range of T* = 30-60 kK. For the latter we also investigate the detailed 
Z-dependence. 

4.1. Hot objects: early-type WC stars 

WC stars are expected to be chemically homogeneous stars in the phase of 
central He-burning. Because of their strong mass loss they show exceptionally 
strong emission lines of He, C and O. They are an ideal target for our model 
calculations because they obey a mass-l uminosity relation which has been de- 
rived theoretically e.g. by lLangerl (jl989). This means that the number of free 



parameters in our calculation is reduced by one. When the luminosity and the 
chemical composition for a specific object are given, e.g. from spectral analy- 
ses, only the stellar core radius R+ or, equivalently, the corresponding effective 
core temperature T* remains. In fact, this parameter is also given by spectral 
analyses. There are, however, striking discrepancies between the observed T* 
(50-90 kK) and the theoretically expected core temperatures from stellar struc- 
ture calculations (120-150 kK). 
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Figure 2. Top: Acceleration within our hydrodynamic WC star model 
in units of the local gravity. The wind acceleration a w i n d due to radiation 
+ gas pressure is in precise agreement with the mechanical + gravitational 
acceleration a moc h + a g rav Bottom: Fe- ionization structure. 



There are in principle two ways how this discrepancy might be resolved. 
First, due to their strong winds, the WC star atmospheres are so extended 
(above the sonic point) that the photosphere is located far out in the wind (at 
~ 3R±). The hydrostatic layers are therefore not directly observable. Second, 
He-burning stars might develop an extended conne ctive envelope below the sonic 
point. Based on static models, llshii et alJ §999) have shown that this might 
be the case for very high luminosities and/or metallicities. Note, however, that 
the extension in these models occurs when the Eddington limit is approached 
close to the 'hot' Fe-opacity peak. Presumably because of the hydrostatic as- 
sumption, these models form an extended convective zone instead of launching 
an optically thick wind. Our hydrodynamic wind models now have the potential 
to discriminate between these two possibilities. On one hand they provide the 
sub-photospheric wind structure, on the other hand they can show under which 
conditions an optically thick wind actually develops. 

The fir st self-consistent PoWR win d models for an early-type WC star are 
described in lGrafener Hamannl (|2005h . These calculati ons are based on a pre - 
vious spectral analysis of the WC 5 prototype WR111 ((Grafener et al.l I20Q2T ) . 
where an effective core temperature of T* = 85 kK was obtained. For the hy- 
drodynamic calculations, the models now had to be extended by the Fe-peak 
opacities, namely the Fe M-shell ions Feix-xvu. First tests however demon- 
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Figure 3. WR 111: Comparison of the synthetic spectrum from our hydro- 
dynamic 13.6 Mq WCE model (dashed line, grey) with observations. 



strated that no wind- driving is possible for 'cool' stellar temperatures around 
85 kK. The models in this temperature regime showed a strong deficiency in the 
radiative acceleration directly above the sonic point. If, on the other hand, T* 
is increased to a value of 140 kK (at fixed luminosity) this problem is resolved. 

Remaining problems in intermediate wind layers (at velocities around 1000 
km/s) are compensated by the choice of a relatively high clumping factor (D = 
50). Wind clumping reduces the observed mass loss rates by a factor of y/D 
because it enhances recombination processes. When the wind density is reduced 
by this factor, the typical WR emission line spectrum remains constant. Our 
models now show that the radiative force behaves similar, i.e., when clumping 
is increased and M decreased the force remains similar which means that a ra d 
increases. Note that iGrafener fc Hamannl (|2003l ) have shown that O star winds 
show the same behavior. 

For the hvdrodvnami c WCE model we assume a stellar mass of 13.6 Mq ac- 
cording to the relation bv lLangerl (|l989l ) for a luminosity of L+ = 1O 5,45 L0. The 
resulting wind acceleration is plotted in Fig. [5] together with the Fe-ionization 
structure. With the obtained mass loss rate of M = lO~ 5 ' 14 M0/yr and the 
terminal wind velocity of Voo = 2010 km/s the observed spectrum of WR111 
is convincingly reproduced (Fig.|3J). Nevertheless, the weakness of the electron 




Figure 4. Same as Fig.|31 but a reduced mass of 12 M Q has been adopted 
in the hydrodynamic calculation. 



scattering wings (e.g. of C iv 5808) indicates that the choice of D = 50 is too 
high. Although the clumping factors can only be determined very roughly, 
WC stars usually show values of the order of D = 10. Our models might there- 
fore underestimate the mass loss rate by a factor of v5, presumably due to the 
incompleteness of opacities. 

Concerning the temperature re gime where the sonic p oint is located, our 
models are in line with the results of lNugis fc Lamersl (|2002| ). However, because 
of the higher effective core temperature in our model, the temperature of T s = 
199 kK is reached at smaller optical depth r s = 5.4 (cf. Eq.EJ). Furthermore, our 
models show an extreme sensitivity to the L/M ratio. A change of the stellar 
mass from 13.6 to 12 Mq increases the mass loss by a factor 1.5. As shown in 
Fig. 01 the resulting model spectrum corresponds to a later spectral type. 



4.2. Cool objects: WNL stars 

Hamann et. al (this volume) have analyzed a large part of the galactic WN 
stars by means of a grid of line-blanketed PoWR models ( Hamann fc Crafenerl 
l2004h . The stars in this sample divide into two distinct groups. The first group 
consists of a mixture of early- to intermediate-type WN stars with luminosities 
below 1O 6 L . The second group consists of late-type WN stars (WN6-9) with 
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Figure 5. WNL mass loss rates for different stellar temperatures and 
masses (indicated by the Eddington factor r c ). The model for WR22 is 
calculated with an enhanced hydrogen abundance (see text). 



luminosities above 1O 6 L0. The latter are located to the right of the main 
sequence (at T+ = 35-55 kK) and show hydrogen at their surface. The extreme 
brightness of these stars already implies high L/M ratios. We have therefore 
calculated a grid of WNL models with a fixed luminosity of 10 6 ' 3 Lq and core 
temperatures in the range of 30-60 kK. For the stellar masses, values of 67 and 
55 M & were adopted, corresponding to Eddington factors T e = XeL-k/'itrcGM* 
of 0.55 and 0.67. Typical WR surface abundances of Xn = 0.2 and = 0.015, 
and a clumping factor of D = 10 are assumed. 

The resulting mass loss rates are plotted in Fig.[5J As expected for optically 
thick winds, the mass loss shows a strong dependence on the ratio L/M or 
equivalently on T e . Moreover, the mass loss increases with decreasing T*, which 
can also be understood in the framework of the optically thick wind theory. 
At cooler T e g- a higher optical depths is needed at the sonic point to reach the 
demanded temperatures (see Eq.|HJ). Therefore higher wind densities are needed. 
The obtained synthetic spectra reflect the observed sequence of WNL spectral 
subtypes, starting with WN6 at 55 kK to WN9 at 31 kK. 

To verify our assumption of a high L/M ratio for WNL stars we have per- 
formed a more detailed investigation of the WN 7 component in the eclipsing 
WR+O system WR22. The Ostar in this system is so faint that the flux distri- 
bution in the optical and the UV is dominated by the WR star. Nevertheless, 
mass estimates are available from orbital solutions for the marginally visible 
Ostar absorption features . On this way values of Mp sin 3 i = 55 db 7 M 
( Schweickhardt et al.lll99flh and 72 ± 3M (|R.auw et al.l fl996) have been deter- 
mined, both with a ratio of M^r/Mo = 2.7. 

With our self-consistent wind models we obtain a reasonable spectral fit for 
WR22 with T e = 0.67, a slightly enhanced hydrogen abundance Xn = 0.4, and 
a stellar temperature of 45 kK (see FigEJ). With the standard distance modu- 




lus of 12.1 mag to Car OBI, our model with L* = 10 63 Lq exactly reproduces 
the observed flux distribution. The corresponding s tellar mass of 78 M r^ is in 
excellent agreement with the mass determination bv lRauw et al. I (Il996h . Note, 
however, that the resulting stellar mass depends on Xh (because T e depends on 
Xh) and on the exact stellar luminosity, i.e. t he adopted distance. Wit h a dis- 
tance modulus of 12.55 mag, as determined bv lMassev &: Johnson! (|1993T ) for the 
Carina OB clusters Tr 14 and Tr 16, we would obtain an even higher luminosity 
and therefore an even higher L/M ratio. In any case, WR22 has a luminosity 
equal or above 10 63 Lq and is located rather close to the Eddington limit. 

4.3. WNL stars: metallicity-dependence 

The WR-type winds in our models are predominantly driven by radiation pres- 
sure on iron line opacities. We therefore expect a strong dependence of WR mass 
loss on Z. To investigate this dependence we have prepared a model grid for 
late-type WN stars which is based on our model for WR22. Because we already 
know about the importance of the L/M ratio, we have varied the Eddington 
factor T e (or equivalently the stellar mass) in addition to the metallicity Z (note 
that we scale all metals including nitrogen with Z). As shown in FigJTJ our 
models indeed show a strong Z-dependence. For models with fixed L/M ratio 
we find rather similar exponents as IVink &: de Koterl (|2005l ) (M cx Z 0,86 ). In 
addition, however, our models show a strong dependence on r e . In particular, 
we find that optically thick winds with high WR-type mass loss rates can even 
be maintained at extremely low Z, if the star manages to come close enough to 
the Eddington limit. 
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Figure 7. Z-dependence for WNL stars: Mass loss rates (top) and terminal 
wind velocities (bottom) as obtained from our grid of hydrodynamic models. 



The obtained terminal wind velocities are shown in Fig(7] (bottom) . First, 
our models clearly predict decreasing wind velocitie s with decreasing Z. This 
findin g is generally confirmed by observations (see e.g. lConti et al Jl98fll : ICrowtherl 
2000), although the observational evidence is not as clear as our prediction. 
Moreover, our models show a rather striking behavior because the terminal 
velocities stay constant, or even increase with increasing mass loss rate. Ob- 
viously, our dense wind models violate th e well-established win d momentum- 
luminosity relation for OB star winds (e.g. iKudritzki et alJll9 99). which would 
imply Mvqq = const. Only at the lowest densities the terminal wind velocity 
starts to increase as expected. A closer inspection of the models shows that the 
changes at high wind densities are related to changes in the ionization structure. 
It therefore seems that the velocity structure is dictated by ionization effects 
rather than by CAK-type wind physics. This is in line with our result that, 
in contrast to standard assumptions, the effective f orce multiplier parameter a 
shows values close to zero in our WC star model (see lGrafener fc Hamannl feoOS). 
Also the present WNL models show a tendency towards low a- values (~ 0.2), 
although not as severe as our WC model. 
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5. Conclusions 

We conclude that models for radiatively driven Wolf-Rayet winds imply a strong 
Z-dependence of the mass loss. In addition, the L/M-ratio turns out to be at 
least equally important. Proximity to the Eddington-limit, or over-luminosity, 
might even be pre-requisite for WR-type mass loss. This would naturally explain 
why He-burning stars and extremely luminous stars tend to show WR spectra. 
Moreover it implies that rotation (which is not yet included in our models) 
might play an important role. First, rotation enhances the L/M ratio in the 
H-burning phase by mixing hydrogen into the stellar core. Second, rotation 
reduces the effective gravity at the stellar surface and thereby also enhances the 
effective Eddington factor. Most importantly, we can state that stellar winds 
from WR stars can now be explained by self-consistent hydrodynamic models. 
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